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month period starting in August 2004. The cluster primary emission is 
from very hot intracluster gas with kT ~ 15 keV, but this component 
does not by itself provide the best fitting model. A secondary emission 
component is quite clearly needed, and while this could also be thermal 
at a temperature significantly lower than kT ~ 15 keV, the best fit (to 
the combined PCA and HEXTE datasets) is obtained with a power law 
secondary spectral component. The deduced parameters of the non- 
thermal (NT) emission imply a significant fractional flux amounting to 
~ 25% of the integrated 3-50 keV emission. NT emission is expected 
given the intense level of radio emission, most prominently from a large 
extended ('halo') central region of the cluster. Interpreting the deduced 
NT emission as Compton scattering of the radio-emitting relativistic 
electrons by the CMB, we estimate the volume- averaged value of the 
magnetic field in the extended radio region to be B ~ 0.4 ± 0.2 fiG. 

Subject headings: Galaxies: clusters: general - - galaxies: clusters: in- 
dividual (A2163) - - galaxies: magnetic fields - - radiation mechanisms: 
non-thermal 
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1. Introduction 

High quality spatially resolved measure- 
ments with current X-ray satellites have clearly 
shown that intracluster (IC) gas is not isother- 
mal. With increasing image detail it will 
likely be found that non-isothermality is a 
general feature of clusters. This is only to 
be expected, given the gas extent, origin, and 
the processes that have affected its evolution. 
The only justification for spectral fitting of 
the emission from a large cluster region by 
a single temperature plasma emission code 
(what has been a standard procedure) was in- 
sufficient spectral and spatial resolution. 

Not only the assumption of isothermality 
is generally invalid when characterizing the 
emission from a large cluster region, the ex- 
pectation that the emission is purely thermal 
may also be doubtful when a wide spectral 
band is considered. This is particularly so in 
clusters with known extended regions of radio 
emission, where the radio emitting relativistic 
electrons give rise to non-thermal (NT) X-ray 
emission by Compton scattering off the Cos- 
mic Microwave Background (CMB) radiation. 
The well known possibility that cluster X- 
ray spectra may have high energy power law 
tails (e.g.,, Rephaeli 1977) has been largely 
ignored, also because of the low sensitivity of 
past X-ray satellites to the detection of the 
predicted low-level NT emission. 

Non- isothermal gas distribution and NT 
emission are of considerable interest both in- 
trinsically, for the detailed understanding of 
the astrophysics of clusters, as well as in the 
use of cluster properties and phenomena as 
cosmological probes. X-ray emission is cur- 
rently our best probe of the properties of IC 
gas; the thermal structure of the gas can yield 
important information on energy exchange 



and transport processes. More globally, the 
detailed gas density and temperature profiles 
are required in analysis of measurements of 
the Sunyaev-Zeldovich (S-Z) effect and its use 
as a cosmological probe of the global param- 
eters of the universe and its large scale struc- 
ture. On the other hand, NT phenomena in 
clusters have the potential to contribute sig- 
nificantly to our understanding of the origin 
of relativistic particles and magnetic fields. It 
is obvious, therefore, that there is strong mo- 
tivation for a more realistic characterization 
of cluster X-ray spectra. 

The search for NT X-ray emission in clus- 
ters has been advanced considerably by the 
RXTE and BeppoSAX satellites. We have 
initiated and analyzed long RXTE observa- 
tions of the Coma cluster (Rephaeli, Gruber 
& Blanco 1999, Rephaeli & Gruber 2002), 
A2319 (Gruber & Rephaeli 2002), and A2256 
(Rephaeli & Gruber 2003). In all three clus- 
ters we found evidence for a second spec- 
tral component. While a secondary emis- 
sion component would generally reflect the 
non-isothermality of IC gas, we found that 
in these three clusters the additional emis- 
sion is likely to be mostly NT in origin. In 
each of the clusters this conclusion is based 
on the relative quality of the statistical fits 
(of thermal vs. NT models), and the high 
physical plausibility of NT X-ray emission in 
the central cluster region in which extended 
radio emission was measured. Similar Bep- 
poSAX searches for NT in these and other 
clusters (including A119, A754, and A2199) 
also yielded evidence for NT emission in at 
least some of the clusters (Fusco-Femiano et 
al. 1999, Kaastra et al. 1999), Fusco-Femiano 
et al. 2000, Fusco-Femiano et al. 2003). Note 
that the significance of the BeppoSAX/PDS 
results has been questioned by Rossetti & 
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Molendi (2003), who claim that the level of 
instrumental error is higher than what was 
previously assumed, and that the detection of 
NT emission in Coma is much less significant 
than originally reported by Fusco-Femiano et 
al. (1999). However, this claim has been dis- 
puted by Fusco-Femiano et al. (2004). 

The moderately distant cluster A2163, z = 
0.203, is one of the hottest, most luminous 
clusters, and a prime target of radio, S-Z, and 
X-ray observations in recent years. Its ex- 
tended central region of radio emission (Her- 
big & Birkinshaw 1994) was extensively mapped 
by Feretti et al. (2001); this radio 'halo' is one 
of the largest and most luminous. Spatial cor- 
relation between bright regions in the Chan- 
dra and radio maps may be interpreted as ev- 
idence for intense merging activity (Marke- 
vitch & Vikhlinin 2001), which is thought 
to enhance the efficiency of particle acceler- 
ation. The gas temperature profile and a 
detailed temperature map were determined 
from XMM (Pratt et al. 2001) and Chandra 
(Govoni et al. 2004) measurements. Here we 
report the results from a very long observa- 
tion of A2163, the most distant of the radio 
('halo') clusters observed with the RXTE. 

2. Observations and Data Reduction 

A2163 Observations with the Proportional 
Counter Array (PCA) and the High Energy 
X-ray Timing Experiment (HEXTE) on RXTE 
were made during 95 separate pointings in 
the period August 24, 2004 - March 11, 2005. 
Data from PCA detectors and 2 were col- 
lected in the 'Good Xenon' spectral mode, 
which produces a 256-channel count spec- 
trum nominally from 2 to 1000 keV. The back- 
ground estimation tool was found to produce 
spectra with a several percent gain error for 
the epoch of these observations; this was cor- 



rected to better than 1% in the software. 
Data from the two independent HEXTE clus- 
ters of detectors were taken in event-by-event 
mode, and were subsequently accumulated 
into 256-channel spectra spanning 12-250 keV. 
To subtract the background, each HEXTE 
cluster was commanded to beamswitch ev- 
ery 32 s between on-source and alternate off- 
source positions 1.5° on either side. 

Standard screening criteria were applied 
to the data segments (Earth elevation angle, 
spacecraft pointing, avoidance of the South 
Atlantic Anomaly, times of geomagnetic ac- 
tivity), resulting in a net exposure time for 
PCA of 528 ks and an average exposure time 
of 212 ks for the HEXTE clusters. The 
HEXTE net observation times were shorter 
than the PCA time because HEXTE spends 
half of the time measuring the background, 
and some more time is lost from electronic 
dead time caused by cosmic rays. 

The PCA background was estimated with 
the 'L7/240' faint source model provided by 
the instrument team, and was corrected for 
errors in detector gain independently in the 
two detectors. Uncertainties in the back- 
ground correction precluded using PCA data 
above 17 keV. The total PCA counting rate 
was 8 count /s (c/s), or 21% of background. 
For HEXTE the net counting rate was 0.6 c/s, 
or 0.4% of background. The PCA background 
level is estimated to be accurate within ~0.5%, 
or 0.2 c/s. The HEXTE background, which 
is directly measured, has been determined 
(MacDonald 2000) to be accurate to within 
a few hundredths c/s in long exposures. 

3. Spectral Analysis 

Inspection of the screened light curves for 
PCA and HEXTE revealed no significant vari- 
ation, as expected for a cluster of galax- 
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ies. Accordingly, we co-added all the selected 
PCA and HEXTE data to form net spectra 
for analysis. Spectra from the individual PCA 
detectors were also coadded, and the small 
differences in gain were accounted for in the 
generation of the energy response matrices, 
following procedures prescribed by the PCA 
analysis team. 

A small energy dependence for systematic 
errors (e.g., Gruber et al 2001) averaging 
about 0.8%, was applied to the PCA spec- 
trum. Additionally, PCA spectral channels 
below 3keV and above 17keV were excluded 
because of sensitivity to artifacts in the back- 
ground model, as well as the rapidly declin- 
ing effective area of the PCA outside these 
bounds. Upon inspection of the high energy 
channels of HEXTE we restricted the analysis 
to data in the range 12-80 keV. 

We fit the joint PCA and HEXTE spectra 
summed over all observations to three sim- 
ple spectral models: a Raymond-Smith (R- 
S) thermal plasma emission model, two R-S 
models at different temperatures, and a R- 
S plus a power law model. In addition we 
have checked whether we can find a meaning- 
ful indication for a temperature structure by 
fitting the data with two R-S thermal com- 
ponents plus a power law. In all cases most 
of the observed flux is in a primary ~ 15 keV 
R-S component. Best-fit parameters and 90% 
confidence intervals are listed in Table 1. The 
best-fit temperature in the single isothermal 
gas model is 15.5±0.9 keV, in good agreement 
with the range determined from XMM (Pratt 
et al. 2001) and Chandra (Govoni et al. 2004). 
The observed ~ 6.7 keV Fe XXV K a line 
yields an abundance of Z = 0.24 ± 0.13 Z Q 
(where Z Q denotes solar abundance), quite 
consistent with previously determined values 
(> 0.2). No cold absorption was measurable, 



and given the 3 keV PCA threshold, none was 
expected. 

The fit to a single isothermal model is only 
marginally acceptable, with \ 2 — 43.2 for 32 
degrees of freedom. Residuals have a high- 
low-high pattern which signals the need for 
another smooth spectral component. When a 
second thermal component is added, best-fit 
parameters are kT\ ~ 20.9 keV, and kT 2 — 
5.7 keV, with the second component account- 
ing for a minor fraction, roughly 1% of the 
3-50 keV flux. For this fit, \ 2 = 35.7, lower 
by 7.5 than the value obtained in the single R- 
S model. The F-test probability of the second 
component is 0.94 for two additional degrees 
of freedom. Although the x 2 is acceptable for 
this fit, the 90% error bounds of [16,80] keV 
for kTi lie outside - in fact, well outside - the 
more tightly constrained values of 11.5 ± 1.5 
keV obtained with ASCA (Markevitch 1996), 
14.6 ± 0.9 keV with ROSAT (Elbaz et al. 
1994), 14.6 ± 0.5 keV with XMM (Pratt et 
al. 2001), and 12.4 ± 0.7 keV measured with 
Chandra (Govoni et al. 2004). We therefore 
view this solution as doubtful. Indeed, the 
range of solutions is very wide because the 
problem is numerically highly degenerate. If 
we consider the possible sets of (kTi, kT2), the 
value x 2 +4.6 defines the joint 90% probability 
contour in parameter space. The lowest tem- 
perature combination permitted in this range 
is kT x = 16.1 keV and kT 2 = 0.8 keV, with 
fractions of the 3-50 keV flux, respectively, 
of 99% and 1%. The highest temperature 
combination has kTi ~ 14.3 keV and kT 2 un- 
bounded, with 13% of the 3-50 keV flux in the 
higher temperature component. These results 
seem to indicate that the two temperature fit 
improves on a single temperature one mostly 
by accounting for emission at the higher en- 
ergy channels. 
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When the second component is a power law 
the best-fit has x 2 — 32.2, fully 11.0 lower 
than for the isothermal fit. The best-fit power 
law index is a = 1.8 with 90% confidence 
bounds of [-2.4, 2.7], and the 3-80 keV flux is 
2.2 • 10 -11 erg/(cm 2 s), with 90% confidence 
bounds of (0.3-4.1) • 10~ n erg/(cm 2 s). The 
F-test probability of this power law compo- 
nent is 0.99 for two additional degrees of free- 
dom. The 3-50 keV flux of the power law 
component is 27% of the total. This com- 
ponent comprises most of the flux at energies 
> 50 keV. The iron abundance is Z = O.32Z , 
with 90% confidence bounds of [0.11, O.79]Z . 
The best-fit fluxes from the combined thermal 
and power law emissions, and that of just the 
power law component, are shown in Figure 1, 
together with the data points. Residuals of 
this best fit model are plotted in the lowest 
panel; for comparison, the residuals to the sin- 
gle temperature model are shown in the mid- 
dle panel. Differences between the latter two 
plots are apparent mostly at high energies. 

The results of the temperature mapping 
with XMM (Pratt et al. 2001) and Chandra 
(Govoni et al. 2004) clearly necessitate in- 
clusion of a temperature distribution in the 
fits. The RXTE FOV includes a region which 
is much larger than that probed by both 
XMM and Chandra. This, and the fact that 
the RXTE lacks spatial resolution, make the 
RXTE insensitive to even large scale features 
in the Chandra temperature map. The above 
fit to a two-temperature model grossly sam- 
ples substantially different temperature com- 
ponents. Just so as to verify that we cannot 
get additional useful information from a more 
detailed spectral modeling, we have consid- 
ered also a model with three components - 
two R-S and a power law. When the model 
parameters are fully unconstrained their val- 



ues are obviously very loosely determined: 
The best-fit to this model - with x 2 = 31.1 
for 28 dof - yields ¥T X = 16.3 ± 1.2 keV, 
a negligible second thermal component with 
kT 2 = 3.7 ± 7.5 keV, Z = (0.28 ± O.16)Z , 
a — 1.3 ± 3.0, and power law normalization 
of (1.3±2.8) • 10~ n erg/(cm 2 s). The decrease 
in x 2 when a power law component is added 
is 4.6, so (from F-test statistics) the signifi- 
cance is about 85% (or, more specifically, 83% 
- 86%, depending on the particular choice of 
'interesting' parameters). However, the \ 2 
hypersurface around the minimum correspod- 
ing to this parameter set is not very deep, 
meaning that a small increase of x 2 above its 
minimum value results in a rather different 
parameter set. In spite of the implied sub- 
stantial uncertainty, we find that the need for 
a power law component is removed only when 
kTi > 30 keV, which we consider to be rather 
unrealistic. 

A power law component at the above de- 
duced level is apparent in the HEXTE mea- 
surements. To check general consistency be- 
tween results from the PCA and HEXTE, 
we have determined the parameters of this 
component by analyzing HEXTE data sep- 
arately. Clearly, due to the relatively high 
(12 keV) lowest energy channel of HEXTE, 
an appreciable fraction of the thermal emis- 
sion cannot be detected by this experiment, 
so we have set the parameters of the isother- 
mal component to its best-fit values from the 
joint (PCA & HEXTE) data. Fitting for 
the power law, we obtain a ~ 1.9lo;g, an d 
F(3-S0keV) ~ (2.5l5;|) x W U erg/(cm 2 s) 
for the power law (photon) index and the in- 
tegrated 3-80 keV flux. Thus, HEXTE data 
alone provide significant additional evidence 
for NT emission whose spectral index has a 
narrower 90% confidence interval than that 
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determined in the joint fits. While this is a 
constrained fit, with only two free parameters, 
it is both physically (by virtue of the antic- 
ipated levels of thermal and NT emissions) 
and practically (given the relative insensitiv- 
ity of HEXTE to emission at low energies) 
well motivated. 

We note that a relatively short ~ 109 ks 
observation of A2163 with BeppoSAX/PDS 
yielded an upper limit on power law emis- 
sion, F(20 - 80keV) < 5.6 x KT 12 erg/(cm 2 
s) (Feretti et al. 2001). Our much deeper ob- 
servation yields a significant detection of NT 
emission with a flux level F(20 — 80 keV) ~ 
l.lijxQ x 10~ n erg/(cm 2 s), whose 90% confi- 
dence interval includes the PDS upper limit. 

4. Discussion 

The results of the above spectral analysis 
clearly indicate that the emission in A2163 
cannot be described by a single temperature 
emission model (as is apparent from positive 
residuals at both low and high energies). A 
better fit is provided by a two-temperature 
model, with the secondary thermal compo- 
nent accounting for the low energy residuals. 
The Chandra temperature map of the central 
cluster region shows a region, roughly 2'X2' in 
area, where the temperature is ~ 9 keV. The 
secondary thermal component we deduce has 
a somewhat lower value, but given the large 
uncertainty in the value we deduce, the differ- 
ence is not significant, as is the rough - due to 
the substantially larger (58' FWHM) RXTE 
FOV, which includes the emission from the 
cooler outer regions - consistency of the rel- 
ative flux values of the two thermal compo- 
nents. 

However, in the best fitting model the sec- 
ond component is not thermal, but rather 
a NT power law. Indeed, NT emission is 
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Fig. 1. — The RXTE (photon) spectrum of 
A2163 and folded Raymond-Smith (kTi ~ 
14.5 keV) and power law (index a = 1.8) 
models. HEXTE data points are marked with 
circles and 68% error bars. The total fitted 
spectrum is shown with a histogram, while 
the lower histogram shows the power law por- 
tion of the best fit. The quality of the fit is 
demonstrated in the lowest panel, which dis- 
plays the observed difference normalized to 
the standard error of the data point. For 
comparison, the residuals of the fit to a single 
temperature model are shown in the middle 
panel. The improved quality of the fit to the 
combined thermal plus power law components 
is apparent mostly at high energies. 
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Table 1: Results of the spectral analysis 



Parameter 


single R-S 


two R-S 


R-S+power law 


X 2 /dof 


43.2/32 


35.6/30 


32.2/30 


kTi (keV) 


15.5 [14.6, 16.4] 


20.9 [15.6, 80] 


14.5[12.1, 21.0] 


Normalization" 


0.450 [0.043, 0.047] 


0.036 [0.001, 0.044] 


0.036 [0.022, 0.046] 


Fe abundance 6 


0.24 [0.12, 0.37] 


0.27 [0.12, 0.54] 


0.32 [0.11, 0.79] 


kT 2 (keV) 




5.7 [0.8, 14.3] 




Normalization" 




0.014 [0.001, 0.045] 


0.0035 


J e (5 keV) (cm- 2 s- r ) 






(1.33±0.45) x 10" 4 


Photon index 






1.8 [-2.4, 2.7] 



Notes: 

All quoted errors are at the 90% confidence level. 
a e.m. = Raymond-Smith emission measure in units of 10~ 14 J N e NndV / 4irD 2 , where D is the luminosity distance 

and N e , Nh are the total number of electrons and protons, respectively. 
b Abundance is expressed relative to solar values. 



expected from the populations of relativistic 
electrons that emit both the extended IC ra- 
dio emission, as well as the emission from the 
several powerful radio sources in the central 
~ 2 Mpc region of A2163. In Table 2 we list 
the measured radio fluxes and spectral (en- 
ergy) indices for the dominant sources in the 
RXTE FOV. In addition to the remarkably 
luminous emission from the most extended, 
regularly-shaped 'halo' region, we include the 
emission from four sources measured by Fer- 
etti et al. (2001). Fluxes and spectral in- 
dices for the 'halo' and 'relic' sources are de- 
termined from VLA measurements at 6 & 20 
cm. To estimate the expected X-ray fluxes 
also from the three tailed (T) sources, we need 
to specify their spectral indices; since these 
were not given by Feretti et al. , we adopt a 
typical value of 0.8 ± 0.1. 

On the relevant spatial scales the radia- 
tion field most tightly coupled to the electrons 
by Compton scattering is the CMB, coupling 
that is only further enhanced in A2163 (z = 



0.203) due to the (1 + z) 4 dependence of the 
CMB energy density. In order to quantify the 
emission that effectively results from the scat- 
tering, we need to determine the spectral den- 
sity distribution of the electrons. Radio emis- 
sion from the central extended source domi- 
nates the overall emission, and since we ex- 
pect the mean, volume-averaged field in this 
extended source to be much lower than field 
strengths in the galactic radio sources, it is 
clear that most of the Compton-produced X- 
ray emission comes from the 'halo'. Nonethe- 
less, we have estimated the emission by rel- 
ativistic electrons in the other sources listed 
in Table 2. 

Compton fluxes of the radio sources in Ta- 
ble 2 can be estimated from their measured 
radio fluxes and estimated mean field values. 
The latter can be determined by assuming 
energy equipartition between particles and 
fields. This assumption may be roughly valid 
in the particles galactic sources, by virtue of 
the relatively short timescales of all the rel- 
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Table 2: Values of radio parameters (and la errors) 



Source 


Flux Density - 20 cm (mJy) 


Spectral (energy) index 


Size (kpc) 


Halo 


155 ± 2 


1.6 ± 0.3 


2070 ± 70 


Relic 


18.7 ± 0.3 


2.1 ± 0.3 


460 ± 40 


J1615-062 (Tl) 


34.5 ± 0.5 


0.8 ± 0.1 


450 


J1615-061 (T2) 


6.0 ± 0.3 


0.8 ± 0.1 


120 


J1616-062 (T3) 


24.9 ± 0.4 


0.8 ± 0.1 


270 



Source size is based on a Hubble constant of 70 km s 1 Mpc 



evant processes governing the particles and 
fields - acceleration, couplings to the field and 
interstellar gas. (We note in passing that the 
attainment of equipartition in the IC space 
is questionable: IC fields and particles are 
likely to be of galactic origin, but the evolu- 
tion during their expulsion from galaxies and 
in IC space is quite different; effective equi- 
libration is not likely to occur there.) Us- 
ing the deduced values of the equipartition in 
the radio sources we have estimated their re- 
spective Compton X-ray fluxes. These turn 
out to be quite small for the relic as well as 
the three tailed sources, whose combined flux 
in the 3-80 keV band is < 1% of the mea- 
sured value. (Note that consideration of the 
dominant non-'halo' radio sources in the clus- 
ters in which we have previously found evi- 
dence for NT emission - Coma, A2256, and 
A2319 - similarly shows that their predicted 
X-ray emission contributes negligibly to the 
deduced NT emission in each of these clus- 
ters.) 

Since the known radio sources in the cen- 
tral region of A2163 do not contribute ap- 
preciable X-ray emission, we attribute all the 
measured emission to electrons in the 'halo', 
and proceed to determine the mean, volume- 
averaged field in this region, B rx , using the 



Compton-synchrotron formulae (e.g., Rephaeli 
1979). For consistency with the basic premise 
that most of the measured Compton flux is 
due to electrons in the 'halo', and the fact 
that the electron spectrum can be more pre- 
cisely deduced from the radio data, we use 
the radio spectral index in order to infer the 
electron index. Doing so we compute B rx ~ 
0.4 ± 0.2 fiG. It should be emphasized that 
the overall uncertainty in the estimated mean 
field is substantially higher than this formal 
error due to the lack of spatial information 
on NT X-ray emission, which necessitates in- 
voking the assumption that the emitting rela- 
tivistic electrons and the field are co-spatial. 

As a consequence of the assumption that 
the spatial factors in the theoretical expres- 
sions for the radio and NT X-ray fluxes are 
roughly equal, it follows that the mean value 
of the deduced magnetic field is independent 
of the source size and distance. The relativis- 
ts electron energy density does depend on 
these parameters. Scaling to the observed ra- 
dius of the diffuse radio emission, and inte- 
grating the electron spectral distribution over 
all energies above 1 GeV, we obtain p e ~ 
1.4±q;| x 10 -13 /ifo erg/cm 3 , where h 70 is the 
value of the Hubble constant in units of 70 
kms _1 Mpc _1 ). Based on the high Galactic 
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proton to electron energy density ratio of cos- 
mic rays, it can be conjectured that the ener- 
getic proton energy density much higher than 
this value. A more specific assessment cannot 
be made without consideration of the electron 
and proton origin and their respective energy 
losses. 

IC magnetic fields can also be estimated 
from Faraday rotation measurements of back- 
ground radio sources seen through clusters, 
yielding a different mean field value, Bf r . 
These measurements usually yield field values 
that £1X6 db few fiG (see, e.g.,, Clarke, Kro- 
nberg, and Bohringer 2001, and the review 
by Carilli & Taylor 2002), up to an order of 
magnitude higher than values of B rx . Aver- 
age field strength of a few /iG over an ex- 
tended cluster region would have important 
implications on the range of electron ener- 
gies that are deduced from radio measure- 
ments, and therefore on the electron (syn- 
chrotron and Compton) energy loss times. 
Weaker mean (volume-averaged) fields imply 
higher electron energies, and therefore shorter 
electron energy loss times, with possibly im- 
portant ramifications for relativistic electron 
models (e.g.,, Rephaeli 1979, Sarazin 1999, 
Ensslin et al. 1999, Brunetti et al. 2001, Pet- 
rosian 2001). However, the apparent discrep- 
ancy between deduced values of B rx and Bf r 
is usually naively interpreted. B rx and Bf r 
actually provide very different measures of 
the mean field strength: B rx is essentially 
a weighted volume-average of the relativistic 
electron density and of B q , with q > 2; on the 
other hand, Bf r is an average of the product 
of the line of sight component of the field and 
the gas density. In general, the overall spa- 
tial dependence of these averages are consid- 
erably different, and so are their deduced val- 
ues. Specific examples and quantitative com- 



parisons, including also the impact of observa- 
tional uncertainties, are given by Goldshmidt 
& Rephaeli (1993), and by Newman, Newman 
& Rephaeli (2002). 

We are grateful to the referee for very help- 
ful critical comments. This work was sup- 
ported by a NASA grant at UCSD. 
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